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By means of free fermionic techniques we study the time evolution of the entanglement entropy, 
S(t), of a block of spins in the random transverse-field Ising chain after a sudden change of the 
parameters of the Hamiltonian. We consider global quenches, when the parameters are modified 
uniformly in space, as well as local quenches, when two disconnected blocks are suddenly joined 
together. For a non-critical final state, the dynamical entanglement entropy is found to approach 
a finite limiting value for both types of quenches. If the quench is performed to the critical state, 
the entropy grows for an infinite block as S(t) ~ In Int. This type of ultraslow increase is explained 
through the strong disorder renormalization group method. 



I. INTRODUCTION 



Recently, we have witnessed increasing interest in 
studying the entanglement properties of quantum many 
body systems 1 - - — in different disciplines: quantum infor- 
mation, condensed matter physics and quantum field the- 
ory. Among the various measures for quantifying en- 
tanglement, the von Neumann entropy and its general- 
izations to Renyi entropies, as well as the entanglement 
spectrum^ have been widely used to obtain useful infor- 
mation about the topological and universal properties of 
an extended quantum system, in particular at a quan- 
tum critical point. For homogeneous, i.e. non-random, 
systems, many basic results are known in one dimen- 
sion from conformal field theory,^ which have been con- 
firmed by exact and numerical calculations on specific 
models.— ~— For a quantum spin chain one generally con- 
siders the entanglement entropy, Si, between a block of 
I contiguous spins and its complement. For periodic 
chains, where the block has two boundary points con- 
nected with the remainder of the system, the entangle- 
ment entropy at the critical point for I 3> 1 scales as 
Si = |ln^, where c is the central charge of the confor- 
mal field theory. Away from the critical point, the en- 
tropy saturates to a value Si = | ln£, where £ <C I (and 
£ ^> 1) is the correlation length of the system. Recently 
universal finite-size corrections to the Renyi entropyi^r— 
as well as the entropy of non single-connected blocks have 
also been studied . 16 ' 17 

If the couplings in the chain are inhomogeneous, 
such as there is an internal defect^ 1 ^— or the in- 
teractions are quasi-periodic or aperiodic,— then the 
prefactor of the critical entanglement entropy, the so 
called effective central charge, c g, is generally dif- 
ferent from that in the homogeneous system. For 
chains with random couplings, c c g has been calcu- 
lated analytically 2 ^— by the strong disorder renormal- 
ization group (SDRG) method^ and numerically by 
frec-fcrmionic technique a 30 ' 31 and by the density-matrix 
renormalization group (DMRG) method^ Also the en- 
tanglement spectrum of random XX chains has been 



studied, both by the SDRG method and numerically^ 
We note that the entanglement entropy can be stud- 
ied even in higher dimensional random quantum systems 
by numerical implementation of the SDRG method^i - — 
provided the critical properties of the systems are con- 
trolled by infinite-disorder fixed points ; 29 ' 37 as in one di- 
mension. 

The nonequilibrium quench dynamics of quantum sys- 
tems has become a very active field of research, both 
experimentally and theoretically^ Dynamical aspects 
of the entanglement entropy are of interest for their 
close relationship to the speed of information propaga- 
tion through an interacting quantum system. In these 
investigations one changes (some) parameters of the 
Hamiltonian suddenly, and asks how the entanglement 
evolves in time^ 9 - One generally distinguishes two types 
of quenches: global and local quenches. For a global 
quench the parameters are changed everywhere in space. 
In this case the entanglement entropy has a linear in- 
crease in time t, irrespective of the initial and the fi- 
nal state of the system. This type of dynamics has 
been explained in terms of quasiparticles (elementary 
excitations)^ In the other type of quench, known as 
a local quench, the parameters of the Hamiltonian are 
changed only locally; for example, a block, which is dis- 
connected from the rest of the system for t < 0, is instan- 
taneously connected at time t — 0. For the local quench 
the entanglement entropy at the critical point is found 
to display a universal logarithmic increase^ Si — ^ In t, 
( C this relation has been later derived through con- 
formal invariancei 41 ' 42 

Concerning dynamical entropy in inhomogeneous sys- 
tems there have been only a few studies in specific 
situations ! 20 ' 32 ' 43 " — For local quench the effect of defects 
has been studied. The defects can be, for example, in 
the form of couplings between a block, where the entan- 
glement entropy is studied, and the rest of system; for 
this case, the prefactor of the logarithmic i-depcndcncc 
of the entanglement entropy is found to be the same 
as measured in the static case.— On the other hand, 
quenched disorder changes the entanglement dynamics 
in a more drastic way. Previously, the entanglement en- 
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tropy dynamics of the disordered Hciscnbcrg chain fol- 
lowing a global quench was numerically studied^ where 
a slow increase of the entropy with time was observed; 
the numerical data in the time regime (i < 500) obtained 
by time dependent DMRG suggested that the entropy 
grows logarithmically with time. The slow propagation 
of signals in disordered system has been later supported 
by theoretical work^ 3 - in which, by means of the gener- 
alized Lieb-Robinson boundj 4 ^ a bound for time evolu- 
tion of the entanglement entropy is derived in the form: 
Se(t) < Ci + c 2 log(£ \t\), with C\ and c 2 being constants. 

In this paper we revisit the problem of the entan- 
glement entropy dynamics in disordered quantum spin 
chains. The model we consider is the random transverse- 
field Ising chain. Our study extends previous investiga- 
tions on entanglement dynamics in disordered systems in 
several respects: (i) we study the entanglement dynamics 
both at the critical point and in the off-critical phases, 
using numerically exact free-fcrmionic techniques; (ii) we 
consider both global and local quenches; (iii) we study 
the time evolution for a very long period of time and 
obtain the long-time asymptotics for finite systems; (iv) 
furthermore, we explain the numerical findings based on 
SDRG. 

The structure of the rest of the paper is as follows. 
In section |TT] we introduce the model, describe its basic 
equilibrium properties and outline the method of calcu- 
lation. Results of the entanglement entropy dynamics 
after global and local quenches are presented in sections 
Illll and llVl respectively. The results are discussed in sec- 
tion El 

II. THE MODEL 

The model we consider is the quantum Ising chain of 
length L defined by the Hamiltonian: 

L L 
i=l i=l 

in terms of the Pauli matrices a^' z at site i. In this pa- 
per we will take periodic boundary conditions so that 
ax+i = <J\. The homogeneous model with the couplings 
J L = 1 and the transverse fields hi = h is in the dis- 
ordered (ordered) phase for h > 1 (h < 1), and the 
quantum critical point is located at h = 142, The crit- 
ical point of the model is described by a conformal field 
theory with a central charge c = 1/2. In the random 
model with quenched disorder, the Ji and the hi are po- 
sition dependent, and are independent random numbers 
taken from uniform distributions in the intervals [0,1] 
and [0, l]h, respectively. The random model is in the dis- 
ordered (ordered) phase for h > 1 (h < 1) and the ran- 
dom quantum critical point is at h = 1. The equilibrium 
critical properties of the random chain has been studied 
by the SDRG method^ and the random quantum criti- 
cal point is found to be controlled by an infinite-disorder 



fixed point, at which the scaling is extremely anisotropic, 
so that the typical length, £, and the typical time, r, is 
related as 

lnr~f, (2) 

with an exponent tf> = 1/2. 

In this work we study the entanglement entropy of 
a block of contiguous spins sitting on sites 1,2, ••■ ,£ 
in the chain; the entanglement entropy is defined as 
S( = Tt( [pi In pi] in terms of the reduced density matrix: 
pi = Tr.;>^|0)(0|, where |0) denotes the ground state of 
the complete system with L sites. In the calculation we 
make use of the fact that the Hamiltonian in Eq. (JTJ) can 
be expressed in terms of free fermions ; 46 ' 48 and the den- 
sity matrix of the corresponding free fcrmionic system is 
then obtained from its correlation matrixj^&i 1 - This cal- 
culation for a system in equilibrium is straightforward. In 
the nonequilibrium case with quench dynamics one has 
two Hamiltonians, say Ho for t < and H for t > 0, both 
in the form of Eq.([T]) but with different parameters. The 
time-evolution of the density matrix is governed by "H, as 
p(t) = cxp(— iTit) p cxp(iT-Lt) , but its matrix elements are 
calculated through the eigenstates of the initial Hamil- 
tonian Hq. Details of the calculation of the dynamical 
entropy can be found in Ref. H^. 

In the following, we first present results for global 
quenches and then for local quenches. In each case, we 
will first briefly discuss the homogeneous model to com- 
pare with our main results for the random chain, which 
will be given subsequently. 

III. GLOBAL QUENCH 

In a global quench the parameters of the Hamiltonian 
are modified everywhere in space. Concerning the Hamil- 
tonian in Eq. ([T]), we modify the transverse fields, but 
leave the couplings unaltered in the quench procedure. 

A. Homogeneous chain 

In the homogeneous chain the transverse fields are 
changed from ho to h and we measure Si (t) in a chain of 
total length L = 256 for various sizes of the block £. Dif- 
ferent combinations of ho and h are considered, including 
quenches from an ordered state to another ordered state 
[Fig.Q](a)], from a disordered state to another disordered 
state [Fig. [1] (d)] and quenches through the critical point 
[Fig. [I] (b) and CD (c)]. 

In each combination of ho and h the dynamical en- 
tropy has a similar behavior. After a linearly increas- 
ing period, S = a(ho,h)t, the entropy saturates to a 
value S = f3(ho,h)£, and decreases subsequently. This 
time dependence of S(t) repeats quasi-periodically, which 
is different from the case in the thermodynamic limit 
L — > oo, where the entropy remains constant after it 
saturates. 39 ' 49 
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FIG. 1: (Color online) Entanglement entropy of a block of I 
contiguous spins in the homogeneous chain of length L = 256 
after a quench from a state with ho to a state with h. Results 
of the semiclassical calculation are given by the gray curves. 



We recall that the exact values of a(ho, h) and /3(ho, h) 
for L — > oo and for a large £ have been calculated^ for 
the special case, when the quench is performed to the 
critical state with h = 1, the coefficient a(ho, h = 1) can 
be calculated in a closed form^S Following a semiclassical 
approach in terms of ballistically moving quasiparticlcs 
formulated in Ref.[5l|, we are able to calculate the dynam- 
ical entropy for the finite chain. The quasiparticles are 
Fourier transforms of kink states, created as a pair of en- 
tangled free fermions with quasimomenta ±p . 39,51 With 
energy e p , a ballistically moving quasiparticle has the 

semiclassical velocity v p = -gjj-. These quasiparticles are 
created homogeneously in space at t = with an occupa- 
tion probability f p . If a pair of entangled particles arrive 
simultaneously in the block and outside the block, they 
will contribute s p = — (1 — f p ) ln(l — f p ) — f p ln/p to the 
entanglement entropy Summing up the contributions 
from all quasiparticles (i.e. over positions and quasimo- 
menta) we obtain the dynamical entropy as shown in 
Fig. [1] by the grey curves. As can be seen, the results 
obtained from this semiclassical approach fit our numer- 
ically exact data perfectly, even for a large range of t. 




FIG. 2: (Color online) Disorder-averaged entanglement en- 
tropy for several global quenches performed outside the criti- 
cal point of random chains, plotted against double logarithmic 
time In Int. (a): quench from a fully ordered initial state with 
ho — > to an ordered state with h = 0.6, 0.8 and 0.9; (b): 
quench from a fully disordered initial state with 1/ho — > to 
a disordered state with 1/h — 0.6, 0.8 and 0.9. System sizes 
ranging from L = 32 to L — 192 are considered; the bigger 
the system, the higher the value of the entropy at large t. 



hi(t < 0)/hi(t > 0) = h /h,\/i. We also consider the case 
where the quench is performed at the critical point, i.e. 
ho = h = 1; for this case we use two independent sets 
of random variables for the transverse fields before and 
after the quench. 

We have calculated the disorder-averaged entangle- 
ment entropy between two blocks of length t = L/2 in a 
chain for different system sizes up to L — 256. We have 
used at least 10,000 disordered realizations to obtain the 
disorder average. 

The time-dependence of the average entanglement en- 
tropy for quenches to a non-critical state and for quenches 
to the critical point is qualitatively different. Below we 
present results for these two cases separately. 



Quench to non- critical states 



B. Random chains 

Now we turn to the disordered chain with random cou- 
plings and random fields. In our quench procedure, the 
set of couplings {Ji} for a given sample remains unal- 
tered, whereas the width of the transverse-field distri- 
bution is changed from ho (t < 0) to h (t > 0). We 
consider quenches in the off-critical phases as well as 
quenches to the critical point. In most of our calcu- 
lations the widths ho and h are different; we keep the 
local transverse field on each site correlated before and 
after the quench, but change the relative magnitude by 



We performed quenches from a fully ordered initial 
state ho — > to ordered states with < h < 1 and 
quenches from a fully disordered initial state ho — > oo 
to disordered states with finite h > 1. Results for the 
disorder-averaged entanglement entropy as a function of 
time for quenches in the ordered phase with parameters 
h = 0.6, h = 0.8 and h = 0.9 are shown in Figl^a). 
Results for quenches in the disordered phase with pa- 
rameters 1/h = 0.6, 1/h = 0.8 and 1/h = 0.9 are 
shown in Figl^b). In both cases the time variation is ex- 
tremely slow, therefore we have used double-logarithmic 
time scales in the figures. 

As seen in the figures, the entanglement entropy for 
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FIG. 3: (Color online) Saturation values of the entanglement 
entropy at t — !> oo and for large L. for quenches from a fully 
ordered initial state ho — ¥ to an ordered final sate, as well as 
from a fully disordered state ho — > oo to another disordered 
state h > 1. The results for quenches in the region h < 1 
are plotted against h, while the results for quenches in h > 1 
are plotted against 1/h. Inset: Finite-size corrections to the 
saturation value of the entropy for a quench from ho — > to 
h = 0.8 and h = 0.9, as well as from l/ho -¥ to 1/h = 0.8 
and 1/h = 0.9. The asymptotic values in the large- L limit 
used for this plot are 5(0,0.8) = 1.615, 5(0,0.9) = 1.941, 
5(c3o,0.8 _1 ) = 1.450 and 5(oo,0.9 _1 ) = 1.823. 



different system sizes L saturates to a value <Sl(/io, h) at 
large t. Furthermore, this saturation value converges to a 
L-indcpendcnt value for large L, i.e. lim£_>. 00 <Sl(/io, h) — 
S(ho, h). This asymptotic value for L — > oo and t — > oo is 
larger if the final state is closer to critical point h = 1, as 
shown in Fig. [3] Furthermore, we have obtained an expo- 
nential relation for the finite-size correction term, given 
by S(ho, h) — Sz,(ho, h) sa exp(— L/£(/i)), as illustrated in 
the inset of FigJ3]for quenches from ho — > to h = 0.8 
and h = 0.9, as well as for "dual process" , starting with 
l/ho -> and ending at 1/h = 0.8 and 1/h = 0.9. From 
this relation we have estimated the values of £(h) for h 
close to the critical point, and obtain £(0.8) « £(1.25) = 
19.8(8), and £(0.9) « f(l.H) = 95.(5), which are in 
agreement with the scaling form of the equilibrium corre- 
lation length of the random transverse-field Ising chain:— 
£(/i) ~ \h— 1| -2 - For this study we have used several ini- 
tial states in the ordered and disordered phases, and the 
same scaling form for £(/i) has been found. 



2. Quench to the critical point 

Results for the time dependence of the average en- 
tanglement entropy after global quenches from a fully 
ordered state (ho —> 0) and a fully disordered state 
(l/ho —> 0) to the critical state are shown in Fig.|4j 

The entanglement entropy for the quench to the critical 
point increases with time up to t(L), after which it sat- 



FIG. 4: (Color online) Average dynamical entropy after a 
quench to the critical point of the random chain, (a): The 
initial state is fully ordered (ho — ¥ 0), (b): the initial state 
is fully disordered (l/ho — s> 0). The dashed lines in both (a) 
and (b) have slope 0.25. 



urates to a value. The asymptotic values of the entropy 
for t > t(L) increase monotonously with L. Analyzing 
the numerical data in Fig. [S]we obtain a logarithmic in- 
dependence: 

S L (h ,l) = s(ho) + b\nL, (3) 

where the prcfactor of the logarithm is found to be inde- 
pendent of the initial state: b w 0.173 « In 2/4; This is 
different than the prefactor for the equilibrium entangle- 
ment entropy: c e ff/3 = In 2/6^ 

For t < t(L), we have found a double-logarithmic 
growth of the entropy in time given by 

S(t) = s + alnlnt, (4) 

with a w 0.25, which is also independent of the initial 
state. This ultraslow growth of the entropy in time re- 
flects the nature of an infinite-disorder fixed point which 
is characterized by the activated dynamics lnr ~ L^ nc . 
Combining Eq. © and Eq. we obtain the exponent 
ipnc = b/a = 0.69(3), which is larger than the exact value, 
tjj = 1/2, known for the equilibrium case. 

We have repeated the calculation for the case when 
the initial state is also critical with ho = 1. A double- 
logarithmic growth of S(t) in time and the same exponent 
ip nc were obtained. 

IV. LOCAL QUENCH 

In the local quench process, we consider the entangle- 
ment entropy of a block corresponding to one half of the 
chain with I = L/2, which is disconnected from the rest 
of the chain, with Jr, = Jl/2 = 5 lor t < 0, and is joined 
up at t = with Jl = Jl/2 7^ 0; the transverse field 
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FIG. 5: (Color online) Saturation value of the average dynam- 
ical entropy as a function of InL, for quenches from a fully 
ordered state ho — ► and a fully disordered state ho — > oo 
to the critical point. The straight dashed lines have the same 
slope In 2/4. 



remains unchanged after the quench. We study the time 
evolution of the entanglement entropy for t > 0. 



A. Homogeneous chain 




FIG. 6: (Color online) Dynamical entanglement entropy of 
the homogeneous system after a local quench versus the 
rescaled time t/L, for L — 256. Upper panel: quench to a 
ordered phase and to the critical point. For the latter case 
the conformal result in Eq. ([5]) is shown by the dashed line. 
Lower panel: quench to the disordered phase. 



interval [0,1]. The two couplings Jl and Jl/2 are re- 
moved for t < 0, and are instantaneously joined to the 
chain at t = 0. For the transverse fields we use the dis- 
tribution described in Sec. II. 



For the homogeneous chain of length L the couplings 
joining the block to the rest of the system for t > are 
Jl = Jl/2 = 1- We have calculated the entropy dy- 
namics after a local quench for different values of the 
transverse field, including the critical value and the val- 
ues for ordered as well as disordered phases; the results 
for L = 256 are presented in Fig. [6] 

At the critical point, h = 1, the entanglement entropy 
oscillates with a period t = L/2, and this periodic func- 
tion can be well described as 
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which was first found in Ref. |20| . and has been derived 
recently through conformal invariance42 

If the quench is performed outside the critical point 
the dynamical entropy grows only up to a finite limit- 
ing value, as shown in Fig. [6j both in the ordered phase 
(upper panel) and in the disordered phase (lower panel) . 
The amplitudes of oscillations of S(t) are reduced for 
large- L and for large t, and the limiting saturation value 
is of the order of 2| ln^om, with the correlation length, 

(hom — |1 — h\~ X i close to the critical point. 



B. Random chains 

For a random chain, the couplings are independent ran- 
dom variables taken from the uniform distribution in the 



1. Quench outside the critical point 

We first discuss the time evolution of the disorder- 
averaged entropy in the non-critical phases. The results 
for two examples in the ordered phase are presented in 
Fig. [7J (a) and [7] (b), and for quenches in the disordered 
phase are in Fig. [JJ (c) and[7J (d). Due to extremely slow 
time evolution, the data for different system sizes L are 
plotted against In Int. 

In the long-time regime, the entropy approaches an in- 
dependent saturation value <Sl(/i), which converges for 
large sizes: liniz^oo <Si(ft.) = S(h). As shown in Fig. [Jje), 
the value S(h) increases monotonously as the critical 
point is approached; furthermore, this asymptotic value 
at each h is smaller than the asymptotic value of the 
entropy after a global quench to an off-critical phase. 



2. Quench at the critical point 

The dynamical entropy after a local quench at the crit- 
ical point is shown in Fig. [8] for different lengths of the 
chain. The overall characteristics of the entropy in this 
figure is similar to that in Fig. |4] obtained after a global 
quench. 

For a fixed length L, there is a characteristic time 
t(L). after which the average entropy is saturated to 
iSl(1). These saturation values follow a logarithmic L- 
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FIG. 7: (Color online) (a)-(d): Entanglement entropy versus 
double logarithmic time after local quenches performed in off- 
critical phases of the random chain for different lengths, (e): 
The data points with filled black circles are the asymptotic 
results for large L and large t in the disordered phase. The 
data in gray are results for global quenches in the disordered 
phase (also shown in Fig. [3]), given for comparison. 



dependence: <Sl(1) = s\ + bi lni, for large sizes [see the 
inset in Fig. [5] . Here the prefactor of the logarithm is es- 
timated as bi w 0.139 « in 2/5, which is smaller than the 
prefactor for a global quench to the criticality: b w In 2/4, 
and is slightly larger than the prefactor of the equilibrium 
entropy: c e ff/3 = In 2/6. 

For t < t(L), the average entropy has a double- 
logarithmic time-dependence: S(t) = si + a\ mini with 
a prefactor a « 0.16(2). Based on the argument for the 
global quench in sec. IIII B 21 we obtain lnr(L) ~ Lfi™ 
with the exponent ^J,° c = 0.87(3); this exponent is larger 
than the exponent i/i nc = 0.69(3), obtained for a global 
quench. 



V. DISCUSSION 

We have studied the time-evolution of the entangle- 
ment entropy in the random transverse-field Ising chain 
after a global and a local quench. The obtained results 
are strikingly different from that calculated for the homo- 
geneous version of the system. The slow dynamics of en- 
tanglement in a disordered quantum chain was observed 
in a previous numerical study^ and has been supported 
by theoretical work.— Our present numerical study pro- 
vides clear evidence showing that this dynamics at the 
critical point is ultraslow and in a double-logarithmic 
form. 
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FIG. 8: (Color online) Disorder-averaged dynamical entropy 
after a local quench at the critical point of the random chain. 
Inset: The asymptotic values at t — ¥ oo (black circles) for 
different system sizes, compared with the static entanglement 
entropy between two halves of the chain (red triangles). The 
black dashed line has slope In 2/5, and the red dashed line for 
the static case has slope In 2/6. 



To explain the difference observed in the time evolu- 
tion of the entanglement entropy in homogeneous and 
in random chains, we apply the semiclassical quasiparti- 
cle picture^ to both cases. In the homogeneous case^i, 
the quasiparticles are related to kinks in the form of 
domain walls between differently aligned spin configura- 
tions, which are created as entangled fermion pairs with 
quasimomenta ±p moving ballistically in the opposite di- 
rections. These entangled quasiparticles will contribute 
to the entanglement entropy between two regions if one 
particle arrives in one of the regions and the other reaches 
simultaneously the other region. From the occupation 
probability of the modes with p, one can calculate the 
entanglement entropy, as discussed in sec. IIII Al and can 
understand the characteristics of the time evolution of 
the entanglement entropy. 

The properties of the quasiparticles and the dynamics 
of the entanglement entropy in the random transverse- 
field Ising chain can be understood from the asymptoti- 
cally exact SDRG. As known from the SDRG, the scal- 
ing properties of the random chain are described by an 
infinite-disorder fixed point, at which disorder fluctua- 
tions are completely dominant while quantum fluctua- 
tions are negligible . 29 ' 47 The ground state of the random 
chain consists of a set of non-overlapping effective spin 
clusters, each of which has a characteristic energy scale 
A c i, given by the excitation energy of the cluster. The 
size of a cluster, t c \, is finite in a non-critical phase, its 
typical value defines the correlation length £. (In the 
ordered phase, there is a giant cluster which is embed- 
ded in finite clusters). At the critical point, where £ 
is divergent, for the largest clusters we have asymptot- 
ically: | In A c i| ~ i\{ 2 [cf. Eq. ©]. The energy-length 
scaling described above is also related to Sinai-diffusion 
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in stochastic processes^— which explains ultraslow dy- 
namics in one-dimensional disordered environments. 

The spins in a cluster defined in SDRG arc maxi- 
mally entangled. Each cluster contributes to the en- 
tanglement entropy between a block and the rest of 
the system by an amount of s c i = In 2 as long as it 
crosses the boundary of the block. In equilibrium, the 
disorder-averaged entanglement entropy at the critical 
point, where the correlation length is divergent, is ob- 
tained by summing up contributions of all clusters in the 
ground state, = ~^2 e <e s c \, yielding Si = ln£ for a 
block of length £^ 

In the noncquilibrium case, the cluster formations in 
the SDRG picture are time-dependent. The time span 
that quantum correlations between different spins to be 
built up is the time for the formation of the cluster con- 
taining these spins, and is given by t c \ ~ A^ 1 . The time 
span td also corresponds to the time in which a signal 
emitted at one end of the cluster arrives at the other end 
of the cluster via a Sinai diffusion. The time-dependent 
entanglement entropy Si(t) at the critical point can be 
obtained by summing over contribution of all entangled 
clusters up to time t: Sg(t) = Ylt <t Scl ~ m £(*) ~ 
Inini, where is a noncquilibrium length-scale, given 
by £(t) ~ ]n{t) 1 '^' n ". In the long-time limit when the 
cluster of size £ c \ < £ are already formed, the entangle- 
ment entropy saturates to a value that is proportional 
to ln£. The explanation through the SDRG description 
holds both for global and local quenches. The main dif- 
ference between global and local quenches is the excess 
energy is finite in a local quench, while it is extensive 
in a global quench. When the excess energy is exten- 
sive, high energy excitations also contribute to the dy- 
namics of entangelement entropy and may be responsible 
for the nonunivcrsal short-time behavior; this is similar 
to the situation in nonrandom gapless systems after a 



global quench: conformal field theory describes the linear 
growth of the entropy, although excited states strongly 
influence the nonequilibrium dynamics^ 

The ultraslow behavior of the dynamical entanglement 
entropy observed in the random transverse-field Ising 
chain should be generic for random quantum systems 
whose critical point is governed by an infinite-disorder 
fixed point; this includes the random XY-chain (which 
has, in fact, an exact mapping with the Ising chain) jil the 
random XXZ-chaiu.^ and the random quantum Potts 
chains Besides, the quantum criticality of the random 
transverse-field Ising model even in higher dimensions is 
also controlled by an infinite-disorder fixed point It 
has been found that at this critical point in higher dimen- 
sions there is a singular contribution to the entanglement 
entropy of the form ASi ~ log £ for a block in a hyper- 
cube form of linear size £^ this singularity is shown to 
be related to the presence of corners. In the dynamical 
process this corner contribution is expected to increase 
in a double-logarithmic fashion. 

Note added: After submitting this paper we noticed 
the preprint by Levine et a&, in which the time- 
dependence of the full counting statistics in a disordered 
fermion system is studied, which is closely related to the 
local quench problem discussed in this work. 
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